The quantum de Finetti representation for the Bayesian Quantum Tomography and 

the Quantum Discord 
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We point out that the quantum de Finetti representation, unique for infinitely extendable ex- 
changeable systems, assigns a non-zero Quantum Discord to uncorrelated systems and thus cannot 
serve as an universal prior distribution in the Bayesian Quantum Tomography. This apparent para- 
dox stems from linearity of the Born rule for the probability assignment in Quantum Mechanics, 
which results in mixing of one's knowledge about the quantum state and the representative of the 
state in one density matrix. 
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Quantum Mechanics is formulated with the concept of 
quantum state, a density matrix p(t) in general, which 
gives the probabilities of outcomes of all possible mea- 
surements on the system. The probabilities are de- 
fined by the Positive Operator Valued Measure (POVM): 
II Q > 0, a = 1,...,K, ^ a II Q = I with the prob- 
ability assignment being linear in the density matrix: 
p a = Tr(nQ,/o), which is the well-known result of Glea- 
son's theorem Q (see also Ref. 0). By measuring the 
quorum of obscrvables one can obtain the set of data suffi- 
cient for complete characterization of the density matrix. 
Feasibility of such a procedure was brilliantly demon- 
strated in the beginning of 1990-es (actually, the 
firstly suggested quantum homodync tomography scheme 
was proved to be such a powerful and efficient tool, that 
the whole field of quantum states/processes reconstruc- 
tion was aptly nicknamed as "Quantum Tomography" 
(QT) 0)- 

Due to irreversible character and statistical nature of 
the quantum measurements, to perform the QT the ex- 
perimentalist needs an ensemble of systems in identically 
prepared states. Then, a statistical estimation procedure 
can be devised allowing one to estimate the density ma- 
trix parameters (e. g. by using the Maximal Likelihood 
Estimation [H, Q or by resorting to the full Bayesian Sta- 
tistical Inference 1 1 0l — 1 1 -4J ] ) . If the number N of systems in 
the ensemble is sufficiently large, the result of estimation 
is expected to converge to the actual density matrix. 

The quantum no-cloning theorem rules out duplicating 
of an unknown quantum state 15| , thus the QT is always 
a process of updating of information about the state (if 
necessary, updating also the parameters of the measure- 
ment device at the same time (l6jh Such an updating 
is the main idea behind the Bayesian Statistical Infer- 
ence method. The convergence property of the likelihood 
function for a generic measurement on the ensemble (i.e. 
the multinomial distribution) to the Dirac delta-function 
in the limit of infinite number of measurements [171 ] (see 



QT was pioneered by K. R. W. Jones [10j, who obtained 
an upper bound on the accessible information obtainable 
from measurement of a pure quantum state, when the 
latter is represented by an invariant prior measure, and 
indicated a measurement scheme (the so-called isotropic 
scheme) which saturates this bound asymptotically. The 
Bayesian approach to the QT for the system consisting 
of 1/2-spins was extensively studied in Ref. [TTJ] , where 
pure as well as mixed quantum states of such spin sys- 
tems were considered. 

We note that measurements in the QT are not re- 
stricted to separate measurements on individual systems 
of the ensemble. For a finite N it is more efficient to 
measure the entire ensemble as a combined system, see 
Refs. [3S]. As single members of the ensemble are 



concerned, it is shown that asymptotically to the order of 
1 /N the effectiveness of the individual measurements ap- 
proaches that of the combined scheme [2l|. However, the 
combined measurement on several systems of the ensem- 
ble allows one to uncover also the correlations between 
the individual ensemble members. 

The central question in the Bayesian QT is the way to 
represent one's incomplete knowledge about an ensemble. 
Using two premises about the ensemble, namely that of 
mutual exchangeability of the individual systems and in- 
finite extcndability of the ensemble, a unique answer to 
this question is the quantum de Finetti representation, 
which follows from the classical de Finetti theorem and 
Gleason's theorem [l2|, 22 1- It also gives the quantum 
Bayesian rule for updating the probability distribution 
[1 31 ] . The quantum de Finetti representation has the fol- 
lowing form 



pW = / d(,(p)p® N , p 



also Ref. [11|]) results in agreement between different 



Bayesian experimentalists. The Bayesian approach in the 



p(g> p(g> p, (1) 
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where dp(p) = dpP(p) and P{p) is the probability den- 
sity. The Bayesian QT then proceeds as follows. The 
experimentalist's prior knowledge about the state is re- 
flected in the prior probability density P(p). The Baycs 
rule for updating the probability density after measure- 
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merits on the first M systems with the data III , . . . , Hm 
(i.e. the POVM measurement results for one or several 
POVMs) reads 

P(p|n 1 ,...,n M ) = — — -— — , (2) 



p(n u ...,n M ) 



where 



p(u 1 ,...,u M ) = J d P p(n 1 ,...,u M \p)P(p). (3) 

Then, the quantum state for the remaining systems of 
the ensemble becomes uM 



P 



(N—M) _ 



d P p(p\n u ...,n M )p 



®N-M 



(4) 



Notice that using the quantum de Finetti representation, 
(HJ) or (U]), one makes predictions also for the correla- 
tions between the individual systems of the ensemble. 
As noted in Ref. [l2j], the exchangeable representation 
(fTj) cannot be carried to the probability theory formu- 
lated in the linear space either over the field of real or 
quaternionic numbers, thus being a unique feature of the 
complex Hilbcrt space. It was also argued that an ex- 
changeable dc Finetti state is a natural substitute for 
the "unknown quantum state" in the QT and that the 
latter has to be banished from it [l2| . 

However, as we show below, the way the quantum de 
Finetti representation accounts for correlations of the in- 
dividual systems of the ensemble bears a considerable 
problem. The problem lies in the fact that the quantum 
de Finetti prior for a quantum state of N exchangeable 
systems almost surely assigns correlations to them, which 
manifest themselves in nonzero Quantum Discord (QD) 
[23|. The QD accounts for non-classical correlations be- 
tween the individual systems of a composite system and 
can be experimentally measured (see for instance, Ref. 
j24l]). But as we discuss below, if it is known that there 
are no such correlations pjj and the information on the 
measured state is limited (e.g. to basic symmetries) 
one cannot combine these two features in a quantum de 
Finetti prior. This problem is even worse: the posterior, 
as given by the exchangeable de Finetti representation 
Eq. @, will also have a nonzero QD almost surely (see 
below). 

Let us recall the QD definition j23|. The QD quan- 
tifies non-classical correlations of two systems A and B 
of a composite system in the Hilbert space Ha ^Hb- 
Given a density matrix p of a composite state, the QD is 
the difference between two versions of the mutual infor- 
mation, D A {p) = I(p) - Qa(p)- One is I(p) = H(p A ) + 
H(pb) — H(p), where H is the von Neumann entropy 
H(p) = — Tr(pln p) and pa,b = ^a,b(p) are the reduced 
density matrices. The other one is defined by optimizing 
over all possible measurements in A and is given as fol- 
lows Qa(p) = H(p B ) -minJ2 k PkH(p B \k), where p B \k = 



Ti A (E k <g> t B p)/Ti(E k ® t B p) is the state of B condi- 
tioned on outcome k in A, and {E^} is the set of POVM 
elements. These two formulations of the mutual infor- 
mation are two quantum generalizations of the classical 
mutual information I (A : B) = (HA) + H{B)-H{A,B). 
On the other hand, the state is of zero QD if and only if 
there exist a von Neumann POVM Hk = IVteK^fcl that 



^(n fc eg) a B )p(n fc ® t B ) 



(5) 



i.e. the p is a state obtainable by a von Neumann mea- 
surement, where only classical correlations remain. The 
QD is currently thought of as a resource for various clas- 
sically intractable tasks including the quantum computa- 
tion [26|428j. For instance, the remote state preparation 
(a variant of the quantum teleportation protocol) based 
on the QD only, i.e. without the quantum entanglement, 
was already implemented experimentally (29j . 

Let us now inspect the QD of the exchangeable state 
([T]) . To this goal we invoke a sufficient condition [3(| for 
non-zero QD of a bi-partite quantum system, which is 
the rank of their correlation matrix R n ^ m - In our case 
the said composite consists of two exchangeable systems, 
i.e. p( 2 ) = £ R n , m A 



Bm , where A„ and B™ are 



bases in the space of Hcrmitian matrices acting in ~Ha,b- 
Then rank(i?) > al A — d B implies D(p) > (note that 
for an exchangeable state the QD is symmetric with re- 
spect to swapping of systems A and B). Consider the 
two-dimensional systems, where the calculations are sim- 
plified with the help of the Bloch vector representation. 
In this case, the unique measure of Eq. (flj can be cast 
as dp(p) = dp(n) = -^dn-ydriidnzPffi), where the Bloch 
vector satisfies n 2 < 1. Hence p(n) = |(1 + na), where 
a = {<7i, o"2) 03} is the vector of Pauli matrices. Then, 
by Eq. (TTJ), the density matrix p^> of system A or B is 



P 



(i) 



dp(ri) — (1 + ria) 



xa 



(6) 



where we have denoted x = J dp(ri)ri. The composite 
density matrix reads 



,(2) - 



1 



p"-' = - ( l<g)t + xa®l + t®xa + ^ 



TijO-i <g> 0~j 



(7) 

with the matrix r defined as r = f dp(ri)ri ® n (in the 
tensor product notation). Now it is a simple observation 
that the correlation matrix R of p^ 2 ' given by Eq. ([7|) 
in the basis (1, a±, 0-2, 0-3) has the following block-matrix 
form 



R = 



1 I x 



ST 



(8) 



hence if rank(r) = 3 > dim('H^) = 2 then D(p< 2 )) > 0. 
But t is full rank for any distribution P(n) provided that 
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it has three-dimensional domain of support in the Bloch 
ball. The simplest example of this class is the point- 
mass distribution with P(n) = Y^aP<x^(™ ~ ^a)i where 
all p a > and e*i, e*2, and being any three linearly 
independent Bloch vectors. 

Next we evaluate the geometric measure of the QD 
proposed for the two-qubit system in Ref. [3(J, which 
gives the distance to the zero QD states in the Bloch 
vector space. It reads 



^max) 



(9) 



where ||r|| 2 = Tr(r T r) (in our case r T = r) and A r 
is the maximal eigenvalue of the matrix A = xx T + ti 
Simple calculations give 



(2)> 



_ min 

4 nif =i 



d/x(ni) / d/z(n 2 )(l + nift 2 ) 



XTll 



i=2,3 



rrij ® rrij 



no 



(10) 



where the vectors rhj, j — 1,2,3, form an orthonormal 
basis in the Bloch space. Eq. ([TO")) shows that almost 
surely D(p( 2 ') > 0, since the matrix in the square brack- 
ets in Eq. (fTOjl is manifestly positive for almost all choices 
of the measure dfi(n), while the scalar factor preceding it 
is also positive. The only exception is the case of a mea- 
sure dp(n) which has support on the one-dimensional 
vector space (parallel to some vector mi) i.e. when P(n) 
is a distribution confined to a line in the Bloch ball. For 
the same reason, the nonzero QD will almost surely pre- 
vail for any finite number of measurements. This is, in 
fact, a quite general feature of the QD, since it was shown 
[3l| that the states of zero QD belong to a zero mea- 
sure subset of all states. Thus, one is led to accept a 
nonzero QD for the result of the Bayesian reconstruction 
for any finite number of measurements, if the de Finctti 
exchangeable density matrix is admitted as a prior. The 
Bayesian experimentalist making measurements on one 
system at a time will not be confused by this, but if a 
more advanced set-up is to be used with joint measure- 
ments on two or more systems at a time, to measure their 
correlations, the problem is bound to arise due to the 
way the exchangeable representation assigns such corre- 
lations. 

The strength of the Bayesian approach lies in selecting 
a judicious prior, reflecting all the information avail- 
able at hand. This is a recurrent theme of the Bayesian 
Statistical Inference in general 32J, l33j . For instance, if 



one knows somehow that there are no quantum corre- 
lations between the individual systems of the ensemble 
(e.g. they are created one at a time), one would like to 
reflect this in the prior information. But since one does 
not know the exact state of the systems in the ensemble, 
e.g. only some symmetry considerations are known for 



the density matrix parameters, one is forced to select a 
straight line in the Bloch ball for the prior distribution 
to have zero QD. But which one should be selected? In 
the current situation the experimentalist has really no 
clues for choosing it. A limited prior information on the 
actual state of the system (a typical case of the QT) and 
zero QD taken jointly do not allow one select a measure 
in representation Eq. ([1]) for the state of the ensemble 
which incorporates all the available information. Thus, 
the exchangeable representation (Qp is not an universal 
prior that fits all (even the most typical) cases. 

It is to be noted that the described problem has rather 
deep roots in the linearity of the Quantum Mechanics it- 
self. The field of Quantum Statistical Inference viewed 
as a variant of the general Parametric Statistical Infer- 
ence introduces one special feature: the linearity of the 
probability assignment on the estimated parameter, i.e. 
the density matrix p describing the quantum state. This 
unique feature forces one to mix one's incomplete knowl- 
edge on the parameter to be estimated and the parameter 
itself. This is seen already on the single system level. 
A prior with the density P{p) for the estimated den- 
sity matrix p invariably leads to a new density matrix 
Pest = J dpP{p)p by the total probability assignment 



d P p( P )p(n a \p) 

Tr(n Q( o est ), 



dpP(p)Tr(IL a p) 



(11) 



where the passage from the first line to the second is pro- 
vided by linearity of Born's rule (and convexity of the set 
of density matrices) with the acceptance of the result p es t 
as the "quantum state" by Gleason's theorem (we note 
in passing that usage of Gleason's result is an impor- 
tant implicit step in the simple and elegant proof of the 
quantum de Finetti representation in Ref. [12|). From 
this point of view, the experimentalist in the process of 
the QT extracts the actual quantum state from such a 
mixture. Indeed, the Bayesian updating in the QT is 
not an unlimited process, but has as a limit the maximal 
possible information obtainable from the ensemble (for 
instance, the Jones limit for the pure state QT with 
the unitarily invariant prior) . After the maximal possible 
information is extracted (to the inevitable imperfections 
of experimental apparatus and restriction to finite num- 
ber of measurements) no further update is possible and 
the experimentalists concludes what is the actual state 
of the ensemble. 

In conclusion, we have shown that the exchangeable, 
i.e. the quantum de Finetti, representation almost surely 
assigns a nonzero Quantum Discord to ensemble of ex- 
changeable systems. If preparation of uncorrelated sys- 
tems is assumed, the simultaneous requirements of zero 
Quantum Discord and exchangeability of the prepared 
ensemble of quantum states do not allow selection of any 
prior at all within the quantum de Finetti representation. 
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Furthermore, we point out that it is the linearity of the 
Born rule for the probability assignments in Quantum 
Mechanics that leads to such contradictory requirements 
since, by Gleason's theorem, one mixes one's knowledge 
about the state and the representative of the state in one 
density matrix. 
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